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Abstract 

Starting from Antonov's discovery that there is no maximum to the 
entropy of a gravitating system of point particles at fixed energy in a 
spherical box if the density contrast between centre and edge exceeds 709, 
we review progress in the understanding of gravitational thermodynamics. 

We pinpoint the error in the proof that all systems have positive spe- 
cific heat and say when it can occur. We discuss the development of the 
thermal runaway in both the gravothermal catastrophe and its inverse. 

The energy range over which microcanonical ensembles have negative 
heat capacity is replaced by a first order phase transition in the corre- 
sponding canonical ensembles. We conjecture that all first order phase 
transitions may be viewed as caused by negative heat capacities of units 
within them. 

We find such units in the theory of ionisation, chemical dissociation 
and in the Van der Waals gas so these concepts are applicable outside the 
realm of stars, star clusters and black holes. 



1 Introduction 

When I first used the concept of Negative Specific HeatQ| (or more correctly 
Negative Heat Capacity) to explain Antonov's remarkable Gravothermal 
Catastrophe || the Statistical Mechanics community thought I was talking non- 
sense. But all astronomers had known since the late 1800s that adding energy 
to a star or a star cluster would make it expand and cool down|3|. Astronomers 
enjoyed this incongruous result, but had not seen it as a real paradox in ther- 
modynamics until ThirringQ emphasised that there is a theorem that specific 
heats are positive. The 'proof as given in Schrodinger's beautiful bookjj| is so 
simple that it is most surprising that it could ever be WRONG. 
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Consider a canonical ensemble in thermal equilibrium 



< E>= Ei exp {-pEk) / ex P (-^) (!) 

i i 

Cv = d -^- = -kf d -^- = kf < <E>? > 

the final expression, which follows from an elementary evaluation of d < E > /d(3 
via (1), is clearly positive, Q.E.D. However, the Astronomers' argument is hardly 
more complicated and involves the kinetic energy T. 

The Virial theorem for a steady state under a potential energy, V, that scales 
like r~ n reads [for external (or edge) pressure p e and volume V = ^Trr^], 

2T + nV = ip e V . 

For gravity n = 1 and the total energy, E, is the sum of kinetic and potential 
parts. So for an isolated gravitational system (p e = 0) 

E = -T < 

but for particles in motion T = |iVfcT so 

dE „ 3 Ar , 

which is clearly negative! 

Luckily ThirringQ saw how to resolve the paradox and I like to think he did 
so in favour of the astronomers but let me give you the historical development 
without first unravelling the paradox. 

Antonov[|| took N particles in a spherical box of radius r e and released them 
with energy E. The total mass was M = Nm. To find the most probable state 
he looked for a local maximum in the Boltzmann Entropy 



S = -k J /ln/ffr 

subject to the constraints 

N — I fd 6 r 



E= / f^—d b T-- / / Gm 2 -^—d b Td' 
J J 2m 2jJ |r-r'| 

The integrals are taken over 6 dimensional phase space, / is the distribu- 
tion function and /' stands for /(r',p'). The density p(r) — J mfd 3 p and 
d 6 r — d 3 rd 3 p, As others had done, Antonov showed that the maximising / was 
a function of e = ip 2 /m — mip where tp is the gravitational potential determined 



2 



self-consistently and this leads to ip obeying the equation for the isothermal gas 
sphere, as tabulated by ChandrasekharQ. We are all used to the entropy being 
very sharply peaked about the most probable state so no-one prior to Antonov 
had bothered to check whether the stationary entropy states were in fact max- 
ima. Antonov found there was indeed a local maximum whenever the radius, 
r e , of the bounding sphere was not too large, and, as r e was increased, the 
density at the edge, p e , dropped compared with that at the centre, p c , due to 
the gravity. However, when p e /p c dropped below 1/709, he found the point 
of stationary entropy ceased to be a local maximum and became a minimax. 
Moreover he was able to demonstrate that there is no global maximum to S at 
any fixed E so that even when maxima exist they are only local maxima. 

Like many original results Antonov's were at first hard to believe. It was 
particularly strange to me that the trouble occurred not when the sphere was 
too small, where the gravity would be large, but when the sphere was too large. 
Nevertheless the effect was due to gravity as it disappeared when gravity was 
absent. To discover what was happening I developed with Wood the thermody- 
namic theory of self-gravitating gas spheres B. In particular we calculated the 
energy E and the heat capacity Cy of isothermal spheres and drew the p e , V 
diagram for adiabats. Figures 1, 2 & 3. 

Figure 1 shows that (—E)r e / \GM 2 ) never exceeds 0.335 for any equilibrium 
gas sphere, whatever the central concentration p c / Pe] so i if a mass M of par- 
ticles is released with energy E (negative) within a sphere of radius r e greater 
than ta — 0.335GM 2 /(— E), there is no equilibrium state for them to go to! 
Furthermore when r e — ta the density contrast is Antonov's 709. 

Figure 2 plots Cy — dE/dT as a function of the density contrast. Cy starts 
at the classical value + |iVfc for a free gas but increases and becomes infinite 
as the density contrast p c /pe approaches 32.2. For larger density contrasts 
Cy starts negatively infinite and increases becoming zero at the Antonov point 
where the density contrast is 709. At this point the entropy ceases to be a 
maximum at constant r e and E, so stability is lost and the dashed prolongation 
of the sequence is unstable (see KatzQ, [| for the proof that they remain so). 
However, isothermal spheres with density contrasts between 32.2 and 709 are 
stable at fixed E and r e and have negative Cy. 

Years of experience with positive Cy systems are not a good preparation for 
this, so let me review some obvious properties of negative Cy systems. 

1. Two negative Cy systems in thermal contact do not attain thermal equi- 
librium - one gets hotter and hotter by losing energy, the other gets for 
ever colder by gaining energy. Thus negative Cy systems can not be 
divided into independent parts each with negative Cy] so negative Cy 
systems are NEVER extensive. 

2. A negative Cy system can not achieve thermal equilibrium with a large 
heat bath. Any fluctuation that, e.g., makes it temporary energy too high 
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will make its temporary temperature too low and the heat flow into it will 
drive it to ever lower temperatures and higher energies. 

3. A negative Cy system can achieve a stable equilibrium in contact with a 
positive Cy system provided that their combined heat capacity is negative. 
To see this imagine that the negative Cy system 'Minus' is initially a little 
hotter (higher T) than the positive Cy system 'Plus'. Then heat will 
flow from Minus to Plus. On losing heat Minus will get hotter (i.e., its 
temperature increases) but on gaining that heat Plus will also get hotter. 
However, because Plus has a lesser \Cy \ its temperature is more responsive 
to heat gain than Minus's is to heat loss. Thus Plus will gain temperature 
faster than Minus and a thermal equilibrium will be attained with both 
Plus and Minus hotter than they were to start with. They also attain 
equilibrium if Minus is initially a little cooler but the reader should think 
that through. 

Notice that this stability is lost as soon as Plus has the same \Cy\ as 
Minus; i.e., when their combined heat capacity reaches zero from below. 

Now we are in a position to explain Antonov's result via a thought experi- 
ment. Imagine a gravitating isothermal gas confined by a sphere of radius just 
less than ta and adiabatically expand the sphere. Work is done by the gas so E 
becomes yet more negative and ta contracts while the sphere's radius expands 
making r e > ta- The inner parts of the gravitating gas are much denser and 
are held in primarily by gravity so the expansion is mainly taken up by the less 
dense gas in the outer parts. Thus the adiabatic fall in temperature of the outer 
parts due to their expansion will initially be greater than the temperature fall of 
the inner parts. Thus there will now be a temperature gradient with the outer 
parts cooler than the central ones. However, as we have seen, the isothermal 
spheres with density contrasts greater than 32.2 have negative specific heats. 
As heat flows down the temperature gradient the central parts contract and get 
hotter while the outer parts held in by the sphere behave like a normal gas so 
they receive the heat and get hotter too. It is now a race; do the outer parts get 
hotter faster on gaining the heat than the inner parts do on losing it? Clearly 
if the outer parts have too great a positive heat capacity they will not respond 
enough and the inner parts will run away to ever higher temperatures losing 
more and more heat to the sluggishly responding outer parts. This is Antonov's 
gravothermal catastrophe. Our criterion for this to happen is that the combined 
heat capacities of the inner and outer parts should reach zero from below and 
Figure 2 shows this happens at precisely Antonov's point. Thus isothermal gas 
spheres with density contrasts greater than 709 arc internally unstable and will 
spontaneously develop temperature differences between centre and edge. Katz 
& Okamoto^) recently emphasised that the effect of fluctuations will decrease 
the 709 limit for simulations without very large N. 

It is of interest to consider what happens at the end of the gravothermal run- 
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away as the centre gets hotter and hotter and denser and denser. If Antonov's 
point particles are replaced by small hard spheres these will eventually get 
so close that they touch. Then the modified system will eventually achieve 
a new higher temperature equilibrium centred on this condensed core. Thus 
we suggested Q that "When the system approaches the gravothermal catastro- 
phe the system undergoes a phase transition in which a core of hard spheres 
in contact with one another is formed". Aaronson & Hansen |l(J later showed 
this was so. The more astronomically important case of a non-relativistically 
degenerate Fermi Gas was likewise predicted to show the same behaviour and 
give white dwarf configurations with total masses below the Chandrasekhar 
Limit |J. For higher mass systems the Fermi Gas becomes relativistically de- 
generate with a pressure-density relationship that is too soft to resist gravity, 
so beyond the Chandrasekhar Limit the Gravothermal Catastrophe leads to 
Black Holes. The fact that these too have negative heat capacities was demon- 
strated by Beckenstein|ll| and confirmed by Hawking's[^2| exact calculation. 
Earlier (1969) in contemplating how quasars lived and died I predicted (Lynden- 
Bell[[l3)) that giant black holes, dead quasars, lay at the centres of most large 
galaxies. Although Maarten Schmidt early|l4| thought this was quite likely, the 
idea gained acceptance only gradually until twenty five years later Miyoshi et 
al.|l5| 1995 found a definitive case in NGC4258. Now Hubble Telescope results 
are widely interpreted as showing giant black holes in many systems. 

These are not the result of the gravothermal catastrophe of stellar dynamics 
(see Section 3). More energy loss via dissipative gas dynamics and radiation 
from quasars is needed to get such large masses into black holes by astrophysical 
processes. 

Thirring's resolution of the paradoxQ we posed in §1 is that two nega- 
tive specific heat systems can not be in thermal equilibrium, so an equilibrium 
canonical ensemble of them is impossible. Thus the supposed 'proof that spe- 
cific heats are positive implicitly assumes the result and in reality only shows 
that extensive s yst ems have positive Cy. 

With Hcrtcl fllq] , Thirring gave an example of a simple system with a negative 
Cy for a range of energies when considered in the microcanonical ensemble. 
When such systems were put in thermal contact the whole showed a positive Cy 
and underwent a phase-transition that corresponded to the region of negative 
Cy in the microcanonical ensemble. 

Some years later (1977) after Van Kampen challenged me over negative 
Cy, my wife and I[|l7j developed an easily calculable gravitational model which 
demonstrated these effects and predicted the temperature of the gravitational 
phase transition in the corresponding canonical ensemble. The model has many 
equal particles confined on a sphere of variable radius r and is governed by the 
Lagrangian 

N 

L = J2 \ mrl (*? + sin2 d <Pi) + \ Mf2 + \GM 2 /r, r <r<r e 

i=l 
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and r is constrained to lie between some small ro and some large r e . The 
particles on the sphere share energy as in a perfect gas. This system has Cy = 
— (TV + \)k whenever a = \GM 2 /(— E) lies between r and r e . However, it 
has positive Cy whenever r is up against one of its limits, i.e., whenever 'a' 
lies outside the range (ro, r e ), which of course happens at very low or very high 
energies, see Figure 4. When we considered many of these whole systems in a 
canonical ensemble the individual spheres were found against one stop or the 
other even when the mean energy was well within the microcanonical region 
of negative Cy. The ensemble now had Cy > and a macroscopic phase 
transition with a latent heat corresponding to the energy difference between the 
spheres on one stop and those on the other. Adding energy (heat) merely led to 
some of the systems formerly against the lower stop, ro, migrating around the 
curve to the higher stop. Thus even macroscopic systems that have negative Cy 
over a wide range of energies behave quite differently when thermally coupled at 
equilibrium. The whole assembly has a positive Cy and a phase transition which 
we determined in detail, see Figure 4. This behaviour led us to ask whether all 
first order phase transitions could be viewed as due to negative specific heat 
systems of molecular size that might cause all phase transitions as discussed in 
the next section. 

However, before leaving this one it is important to discuss what happens in 
practice as well as what happens 'at equilibrium' especially as the metastable 
superheated or supercooled regions in our model can be large. These metastable 
regions have the sphere against one stop or the other. Even in the canonical 
ensemble it takes a gigantic fluctuation to take such a metastable system over 
the hump. It has to gain a thermal energy comparable to the mean energy of 
each whole system, i.e., it needs a fluctuation -y/TV times the typical one. As 
emphasised by Parentani, Katz & Okamoto|ll| for the black hole case, such 
fluctuations will occur almost never - i.e., once in exp(— TV) independent trials. 
Thus to get the thermodynamic phase transition TV can not be large and even 
for N = 10 one must wait for ~ e 10 complexions of the system. Thus these 
metastable regions of our large systems will be exceedingly stable and no phase 
transition of the canonical ensemble will be observed until the system nears the 
top of Figure 4a. However, once on the unstable branch the system will evolve 
along the downslope on the timescale of thermal diffusion and will then become 
so much cooler than the ensemble that on a similar timescale it climbs the other 
branch to regain the 'tops' temperature of (TV + \)kT ~ \xo- If the ensemble is 
cooled, likewise the transition will not be observed until the temperature reaches 
the much lower temperature of (TV + \)kT ~ \xe when the transition will again 
go violently. Thus there will in practice be an enormous hysteresis compared to 
which 'boiling with bumping' will look quite petty. 

Kiessling[[l9] has emphasised that statistical mechanics done properly gives 
transitions like AC rather than ones associated with the ends of the metastable 
regions but for transitions in these macroscopic systems the path AC would 
not be encountered in practice nor normally in simulations unless the numbers 
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of particles involved are small. However, the metastable regions in microscopic 
systems can be surmounted by fluctuations with smaller energy changes so ice 
melts at the temperature at which water freezes rather than at a significantly 
higher one. 



2 A Different View of the Cause of Phase Tran- 
sitions 

In the last section we saw how a region of Cy < in individual systems gave rise 
to a phase transition when many such systems were placed in thermal contact. 

In reference 13 we asked whether all phase transitions (or at least all first 
order ones) can be viewed as caused by negative specific heat elements at a 
molecular level. Can one look at systems that undergo changes of state and 
identify such negative-heat-capacity elements? We shall demonstrate their ex- 
istence for simple models of chemical dissociation or ionisation equilibrium and 
in systems like the Van der Waals gas. 

We consider two types of particles, A and B which are in the ground state 
tied together in pairs AB. We take the forces between them to saturate when 
such a pair exists and shall model the interaction potential to be a 8 function 
so a single bound state exists with binding energy x- Let uab be the number 
density of pairs in some initial ground state. We shall model this system by a 
microcanonical ensemble of boxes each of volume L 3 = n^ B and each containing 
one particle of type A and one of type B. We shall show that each of these 
microsystems shows a negative specific heat when the energy of each box is just 
greater than the dissociation energy. Taking that as our zero point for energy 
we find the number of energy levels < E for a system of two independent free 
particles is 

f E 3 L 3 E>0 
{ ' E < ' 



Tf(E) 



oc 



while for a bound pair the number of energy levels is Tt,(E) oc \(E + x)L 2 ] 3 ^ 2 . 

The total number of energy levels < E is thus t{E) = t/ + t&. Gibbs 
gives both Si — fcln (dr/dE) and 5*2 = fcln(r) as possible expressions for the 
entropy for small systems. We earlier jlTj showed Si gave negative specific heat 
just above the dissociation point so here we consider S*2. Evidently 

kT = k _dE t (E + xf' 2 + AE 3 



dS 2 t' l[(E + xy/ 2 + 2AE2] ' 

where A oc L 3 ! 2 . It is not hard to show that 

k _ dkT / t 
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is negative for some (large enough) values of L for energies greater than zero and 
on up to about 0.8% for some L. Thus if the system is in a large enough box, 
i.e., rare enough, too many dissociate as the energy is increased so the mean 
kinetic energy per free particle decreases because of the energy soaked up by the 
dissociation. In this sense negative Cy elements are associated with chemical 
dissociation reactions. Much the same arguments hold for ionisation reactions 
at higher temperatures. 

Perhaps of greater interest is the demonstration of similar phenomena in the 
Van der Waals gas and the following demonstration is due to 
R.M. Lyndcn-Bcll my wife. We are all familiar with the dip and hump in the 
isotherms of the Van der Waals gas and how Maxwell's construction replaces 
them with the constant pressure phase-transition that occurs in practice. This 
is what will happen in an extensive system with very many identical subsystems. 
However, at or close to the molecular level the extensivity breaks down because 
such tiny systems can not be readily divided into two equivalent pieces. The 
Van der Waals equation includes such effects in its molecular volume term b and 
its mean molecular attraction term which reduces the effective pressure by ap 2 . 
We shall assume that at a molecular level the tiny indivisible elements of Van 
der Waals gas actually obey Van der Waals 's equation and it is only a cooper- 
ative effect of many of them in a canonical ensemble that makes the ensemble 
obey the Maxwell construction and give the phase transition. We saw that this 
is precisely what happened when we took our gravitating systems in Section 1 
and put many of them into a canonical ensemble at equilibrium. What we now 
show is that an element that obeys the full curve of the Van der Waals equation 
inevitably has a negative C p . The hump in the Van der Waals isotherms at 
temperatures below the critical one gives two stable states AB which share the 
same temperature and pressure (the third is unstable). Thus 



For this to happen (dp/dV)r can not have the same sign everywhere so we need 
an unusual region with it positive. 



so normally, with C p positive, InT will increase with S at constant p. However, 
the states A and B share a common isotherm and are at the same pressure so 



so C p 1 = (dlnT/ dS) p must reverse sign and the graph of InT against S can not 
be monotonic. We have, therefore, demonstrated that for any system which has 




Now 
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a hump in its isotherm like that in the Van der Waals gas there is a region of 
negative C p . This is our prime result! 

One learns at school that for perfect gases C p — Cy = R so one might expect 
that systems with C p negative would also have Cy negative. This is not so, 

Cp ~ CV=T {^) p - T {^)y= T {w) T {^f 

dp\ (dV\ = fdp\ fdV 
dT) v \dT) p \dV) T \dT 

Thus C p is only greater than Cy at positive T when, as is usually the case, 

(^§y^j is negative. But when the pressure along an isotherm increases with 

V, as in a phase transition region, C p is actually less than Cy. It is, therefore, 
a moot point whether Cy has to be negative. Here Van der Waals's gas is of 
great interest as an example because for it gives Cy = fiVfc which is positive 
everywhere! Thus the hump and trough in the isotherm is associated with 
negative C p but not with negative Cy. 



3 Aftermath of the Catastrophe 

We earlier left our gas just undergoing the gravothermal catastrophe with the 
central part contracting and getting hotter while the outer part could not raise 
its temperature fast enough to keep up. Since the centre is now denser the 
Antonov point at which the density reaches p c /709 is now inside the system 
rather than at its edge. Thus this point moves inwards through the mass. In 
astronomy we are primarily interested in applying this theory to a 'gas' in which 
each molecule is replaced by a star so that we deal with a star cluster or the 
nucleus of a galaxy. In such systems the timescale for exchanging the 'heat' 
of random stellar motion is shortest at the centre and can become very long 
in the outer parts. As the Antonov point moves in, the timescale of heat flow 
becomes shorter and the regions beyond the Antonov point are too sluggish to 
respond other than adiabatically. The gravothermal catastrophe occurs again 
and again nay continually, with ever higher densities and temperatures at ever 
smaller scales ]2(], gl], |2^] and the precise form of the initial conditions is rapidly 
forgotten. Since the same process is occurring at ever smaller scales we expect 
a similarity solution with the density of the form 

p(r,t) = p c (t)p*(r*) 

where r* = r/r c (t) and r c (t) the core radius will be some fixed fraction of the 
Antonov radius ta where the density is p c /709. Since the halo is so sluggish 
that it is left behind in the every quickening evolution of the centre we can put 
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dp/dt = at large r. Hence p c p* — —r+pcp'^ = and so we may separate the t 
and dependences to get 



r*p* r cPc 

= = —a . (large r) 

P* r c p c 

Since the * variables are independent of the t variables a must be a constant. 
Thus = Ar~ a at large r* and p c oc r~ Q for all t. 

Now the evolution of the system is due to the heat transport whose rate 
is determined by oc y- where T r is the relaxation time within the Antonov 
radius; for it the standard formula is T^T 1 oc (87rG 2 mln N) where v 2 is the 
velocity dispersion of the stars. In our self-similar collapse 

v 2 oc Gjw2/r B « G Pc rl oc Gp\~ 2 l a . 



Hence 



pc i-i 

OC pc - 

Pc 



On integration we find|2l], |23| 

p c oc \t Q -t\- 2a ^ 6 ^ so r c oc [to - t| 2 /( 6 -«) . 

Thus the core radius becomes zero and the central density formally — ► oo at to. 
The core mass is 

M c oc Pc r\ oc |t - t |3(3-a)/(a-a) 

and the core velocity dispersion behaves as 

v\ oc GM c /r c oc \t -t\^ 2a ^^ . 

We must still determine a. 

Now at large distances oc r~ a and the asymptotic form of the constant 
temperature (isothermal) sphere is p oc r~ 2 so for an outward temperature de- 
crease a > 2. However, if p oc r -2 5 there is an infinite binding energy near 
the centre and to create this in finite time would need an impossibly large 
heat flux. Thus from these general arguments 2 < a < 2.5. Detailed calcu- 
lations give a as an eigenvalue and for stellar dynamics all investigators get 
a = 2.22 ± .01 which gives the rather weak dependence v c oc \to — t| - ( a-2 )/( 6 ~ Q ). 
From these one finds M c oc v c 2( ' 3 2 ^ where the exponent is 7.1 for a = 

2.22. If in core collapse v c increased from 300 km/s to 300,000 km/s correspond- 
ing to a black hole then M c would decrease by a factor 2 x 10 21 leaving much 
less than one star so these conditions are not attained. In fact if M c decreased 
by 10 71 then v c would increase by a factor 10 only. However, in fact a new 
phenomenon sets in a high star densities that produces a delightful new twist in 
the evolution. Henon|24J in his most percipient early work suggested that the 
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formation of binary stars would occur at very high densities and that this would 
produce a new energy source. Heggie]25|] did the seminal work on the binary 
creation rate and Betteweiser & Sugimoto[^6| put this into their code for core 
collapse. As Sugimoto had surmised [^7|, binary formation sets in quite suddenly 
during core collapse and releases energy near the centre. Since the core is of neg- 
ative heat capacity it immediately expands and becomes of lower temperature 
than its surroundings. The stage is now set for the inverse gravothermal catas- 
trophe. As heat is conducted into the core it grows and gets colder. The process 
goes on until eventually the core grows so large that it contacts the cooler outer 
parts of the halo. Then the whole core becomes isothermal and eventually feels 
the heat loss to the outside so the gravothermal collapse stars again. A number 
of these giant thermal pulses occur in simulations of post-collapse core evolution 
so it is hard to tell whether the self-similar evolution |23| predicted on the basis 
of a singular core from which a flux of energy emerges due to continued energy 
emission from binaries gives even a roughly correct average evolution. 

Padmanabhan|2^| has given a longer more quantitative review of gravita- 
tional thermodynamics which discussed many of the phenomena considered here. 

Meylan & Heggie[^9] have an up to date review of the more astronomical 
aspects. 

Miller & Youngkins|30| |3lJ have given the mean field theory for a simple sys- 
tem of spheres and have conducted numerical simulations to verify the existence 
of gravitational phase transitions. 

Both Miller & collaborators |3^] and Tsuchiya|33) have conducted very long 
term integrations to investigate the approach to equilibrium of a one-dimensional 
system of mass sheets. For these relaxation is particularly slow as the interac- 
tions when they cross are rather smooth. 
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Figure 1: The dimensionless binding energy -~Er e /(GM 2 ) of an isothermal 
gravitating sphere of mass M in a spherical container of radius r e plotted as 
a function of the density contrast, lnp c / p e = — i>i, between centre and edge. 
Instability sets in at A, the maximum. 



14 




Figure 2: The specific heat Cv = dE/dT for a self- gravitating sphere in a 
spherical container plotted as a function of the density contrast. Instability sets 
in at A as C v reaches zero FROM BELOW! 
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Figure 3: The adiabat plots the surface pressure of the sphere against its volume 
V = ^Trr^ at constant entropy, instability sets in at A for an isolated system 
but at B if the system is held at constant pressure. 
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Figure 4: (a) T as a function of E for the variable radius slippery sphere model. 
AC is the phase transition predicted for a canonical ensemble but replaced in 
practice, if ./V is large, by lines to the maximum (T increasing) or minimum (T 
decreasing) of the T(E) curve which then shows a hysteresis loop, (b) (fcT) -1 
as a function of E. The shaded regions are of equal area for the canonical phase 
transition, (c) The mean radius of the sphere as a function of energy. 
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